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Abstract: We compute the index and associated spectral density for fluctuation operators 
which are defined via the Lagrangian of N = 2 SQCD in the background of non-abelian 
confined multimonopoles. To this end we generalize the standard index calculations of 
Callias and Weinberg to the case of asymptotically nontrivial backgrounds. The resulting 
index is determined by topological charges. We conjecture that this index counts one 
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1. Introduction 

From their discovery on [1, 2], non-abelian confined monopoles and vortices received con- 
siderable attention. The reason for this is that these objects may allow the description 
of non-abelian confinement as an electric-magnetic dual Meissner effect, as it was envi- 
sioned by 't Hooft [3] and Mandelstam [4]. We refer to [5, 6, 7] for recent reviews of these 
developments. 

In [8], among other things, the full perturbative quantum energies and a central charge 
anomaly for confined multimonopoles in J\f = 2 SQCD have been computed. Central 
to this quantum computation is the spectral density of certain operators, which can be 
obtained from an index theorem. The present investigation is devoted to the derivation 
of these quantities. The operators in question are the fluctuation operators obtained from 
the expansion of the Lagrangian around the classical background fields that satisfy the ^ 
BPS equations which describe confined monopoles (and many other field configurations). 
These operators, describing the fluctuations of the full second order equations of motion, 
differ from the fluctuation operators that describe the fluctuations of the first order | BPS 
equations. We will discuss this in detail. 

The index has been computed for many different topologically non-trivial backgrounds 
in different models, in the given context for example for vortices and domain walls [9, 
10], but not for confined monopoles. The difference in the case of confined monopoles is 
that the usual techniques developed in [11, 12, 13] cannot be applied directly. Confined 
multimonopoles depend on all three spatial coordinates, but the nontrivial field dependence 
is essentially concentrated in flux tubes, i.e. the vortices that confine the monopole and 
emanate from it, see figure 1. Consequently, the background fields of the fluctuation 
operators do not fall off asymptotically, or terminate in a vacuum configuration. We 
describe in the following sections how to resolve this problem and develop a strategy that 
generalizes the established methods [11, 12, 13] of the (open space) Callias index theorem 
to the case of asymptotically nontrivial backgrounds. 

Following [11], the index of such operators can be written as the sum of an anomaly 
and a boundary or surface term. Both contributions can be conveniently computed as 
appropriate limits, however, for this to be true for the surface term the background fields 
have to be asymptotically trivial. The latter does not apply to the field configurations 
considered in this paper. For the concrete case of | BPS confined monopoles, the principle 
form of the field configurations is depicted in figure 1. The given geometrical setting defines 
the surface term on a boundary of the form of a cylinder at infinity. Hereby, one has 
especially at the discs D± highly nontrivial field configurations, notably multiple vortices 
which are not even known analytically. We will show how to reformulate these surface 
contributions in the form of an index, though of a generalized form. Its computation can 
again be reduced to an anomaly and a surface term on the boundary dD± of the discs. 
It turns out that only the anomaly on the discs gives a non-vanishing contribution, which 
however depends on the IR-regulator and thus leads to a non-vanishing spectral density. 
Based on some explicit examples we conjecture that the resulting index counts one quarter 
of the real dimension of the moduli space for general confined multimonopoles. 
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Figure 1: The left figure shows a confined monopole in the close up with characteristic scales 
for the monopole size, displaying approximate Coulomb monopole behavior, and the size of the 
confining vortices. The right figure depicts the spatial boundary, a cylinder at infinity, where the 
asymptotic behavior of the (multiple) confined monopoles is specified (see main text) . The cylinder 
at infinity consists of the two discs D± at x 3 — > ±00 and the cylinder wall 3oo at r — > 00, with r 
being the cylindrical radial coordinate. 

The paper is organized as follows: In section 2 we introduce the necessary details of the 
model, the BPS equations and the associated fluctuation operators. In section 3 we briefly 
review the conventional techniques, introduce the necessary generalizations and compute 
the index as a function of the background fields. We also relate the resulting index to 
topological charges and formulate a conjecture regarding its relation to the dimension of 
the moduli space of confined multimonopoles. Finally we discuss the relation between the 
fluctuation operators obtained from the Lagrangian and those obtained directly from the 
BPS equations. In section 4 we summarize our results and give conclusions. In particular, 
we give an outline for the general recipe to compute the index in the presence of asymp- 
totically nontrivial backgrounds. In the appendices we give several details for the proofs 
of the statements in the main text. 

2. Confined Monopole BPS Equations 
2.1 Model and vacua 

The full Lagrangian for M = 2 SQCD with generalized Fl-terms and the conventions that 
we use are given in [8]. Here we need only the bosonic part of the Lagrangian which is 
given by 1 

A,os = - Jr Tr { \F% + \D^ + §[0, <ffl + \ } 

- \D^S lT \ 2 - S({<fl +mi,<j> + mijSu . (2.1) 

The complex scalar field eft transforms in the adjoint representation of the gauge group 
U(N), whereas the Nf SU(2)r doublets Su are scalars that transform in the fundamental 

1 Fot the metric we use the east coast convention n^y = ( — ,+,+,+) and generally summation over 
repeated flavor indices i — 1, . . . , Nf is implied. The covariant derivative is defined as D M = 9 M — iA^ for 
fields in the representation R. The field strength is given by _F M „ = d^A v — d v A^ — A v \. The positions 
of SU(2)r indices I are changed by complex conjugation. 
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representation of U(N) and the flavor symmetry SU(Nt). The latter one is explicitly 
broken to a subgroup by the masses. The bracket in the last term indicates the anti- 
commutator of the given matrices. We are interested in static \ BPS solutions and therefore 
we choose the masses to be real (and ordered): rhi = mj and m\ < ... < mjsr f . See 
[8] for a detailed discussion. The auxiliary field D is an SU(2)r vector which is in the 
adjoint representation of the gauge group. We will keep the auxiliary field for notational 
convenience but it is understood to be on-shell, i.e. 

D = ^(T I J S iJ ®Sf-£t) , (2.2) 

where the tensor product notation is defined w.r.t. the gauge group structure carried by 
the fields. 

The SU(2)r triplet £ of Fl-terms explicitly breaks the SU{2)r symmetry and can 
therefore be rotated to a convenient choice [7]. For the rest of the article we will assume 
that it points in the positive three-direction which we parametrize as 

f= (0,0,1^)*, (2.3) 

with v being a positive constant. 

For the generators of the gauge group U(N) we have the following conventions: The 
hermitian generators are {T A } = {T°,T a }, where {T a } forms ansu(iV) algebra, and satisfy 

[T a ,T b ] = if abc T c , T° = -^1=1 , Tr{T A T B } = \5 AB , (2.4) 
where f abc are the real and totally antisymmetric su(iV) structure constants. 

Color-Flavor locked vacua. The Lagrangian (2.1) provides a particular set of vacua 
which preserve a diagonal subgroup of the gauge group U(N) times the unbroken flavor 
group Hp C SU(N) ^-t- SU(Nf). 2 Up to gauge transformations the vacuum is specified by 
the following vacuum values of the scalars: 

0o = -diag(mi,...,m JV ) , [STTi = vS n i, (2.5) 

where the fundamental scalar SJ 3,0 is written asaJVx Nf matrix, whose entries are zero 
for i > N. These vacuum scalars preserve certain subgroups of the original symmetry, i.e. 

H C 0o H^ 1 = cf>o , H c+F ST HcIf = Sl ac ■ (2-6) 

The transformation Hc+f acts from the left as a global U{N) gauge transformation, and 
from the right as an Hp flavor transformation. Hence, one has Hc+f C d\a,g{U{N) x 
SU(Nf) ) and such vacua are called the color-flavor locked phase. 
The associated symmetry breaking pattern is given by 

U(N) x SU(N f ) ^ U c (l) xH c xH F ^ H c+F , (2.7) 

2 The focus is here on the part of the unbroken flavor group that acts nontrivial on the to be considered 
background fields which carry flavor index i — 1, . . . , N . 
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where we assumed that all masses are of the same scale and that m 3> v. If the first N 
masses form q groups of n r degenerate masses the surviving symmetry group is given by 
[14], 

H c+F = S(U(n 1 )x...xU(n q )), (2.8) 

with ^ r n r = N. It supports monopoles with typical size 1 /Am, the inverse mass differ- 
ence, that are confined by flux tubes of width l/gv, see figure 1. 

2.2 J BPS equations and asymptotics 

We denote the classical background fields that satisfy the BPS equations given below by 

Af e = A k , < f >hkg = <p = ^, S^ = Ei, D 3 hkg = D, (2.9) 

where the index k = 1,2,3 runs over the three spatial directions. The covariant derivative 
w.r.t. the gauge field A k will be denoted by T> k . All other fields are assumed to be zero 
for the classical solutions of interest, in particular one has Sj>at = 0. This ansatz implies 
also that classically the adjoint scalar eft is hermitian. 

With the particular choice (2.3) for £, the first description of confined monopoles in 
the given context and the associated novel BPS equations were given in [2], derived using 
the Bogomolnyi trick. It was shown that the energy density is (locally) minimized if the 
following first order equations are satisfied: 

B k -V2V k <p + 5 k3 D = with: D = ^ ( Sj <g) % - v 2 ) 

V z Vi = 0, V w Ei = 0, (2.10) 

where we have introduced complex coordinates such that T^z — ^i ~t~ 2X^2 and T)w — 
P3 + iT>4, with P4 = i 1/2(9? + m i) when acting on the i'th fundamental scalar Ej (in the 
adjoint representation the mass rrii drops out). The chromo-magnetic field is defined as 

These equations a priori look overdetermined, but as was noted in [15] the equations 
for k = 1,2 are identical to the integrability condition of the last one. This novel set 
of equations offers a plethora of non-trivial field configurations carrying various topolog- 
ical charges [10]. For us the main focus lies with configurations that describe confined 
monopoles. In [16] an approximate solution for a single confined monopole was given for 
the case of £7(1) x SU(2) gauge group and Nf = 2 (though for a different choice of £ ). 
The relation of these equations to the supersymmetry algebra and the associated (tenso- 
rial) central charges were discussed in detail in [8]. Equations with less supersymmetry for 
intersecting vortices were derived in [17]. 

Classical Asymptotics. We consider those solutions of (2.10) that describe multiple 
confined monopoles of the form as depicted in figure 1, such that the field configurations 
have an axial orientation in the x 3 -direction. The main input for the index calculation will 
be the asymptotic behavior, i.e. the topology, of these solutions. 

The axial orientation of the field configurations implies that the asymptotic boundary 
has the form of an infinite cylinder, see figure 1, and accordingly we have to specify the 
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asymptotic behavior: i.) At x 3 — > ±00 the boundary is given by the infinite discs D± at 
which the fields behave like multiple vortices, though in general different vortices at £)+ and 
ii.) For r — > 00, with r being the radial cylindrical coordinate, the boundary is given 
by the cylinder wall 3oo at infinity. The flux is confined in the vortices which are infinitely 
far away from the cylinder wall 3oo and therefore vanishes exponentially with correlation 
length proportional to gv, see (2.10). Hence, at the cylinder wall one has asymptotic vortex 
behavior, with winding in the Higgs fields and the long-ranged gauge fields. Due to the 
monopoles this winding depends on x 3 , but this dependence is exponentially located at 
the monopoles with the characteristic length given by the associated mass difference |Am|. 
Concretely, the asymptotic field behavior is as follows: 

Cylinder Wall 3oc>: The confinement of the monopoles/flux by vortices implies Bk « « A r , 
where ~ means equal up to exponentially suppressed terms ~ e~ 9vr . The Higgs fields Sj 
approach their vacuum values (up to winding) exponentially fast, i.e. with the spatial 
angular coordinate 9 one has, 

» USlf c , U = e ie{w+wc + F ^ . (2.11) 

The part of U which lies in the unbroken color-flavor symmetry Hq+f (2-8) is generated by 
wc+f{x 3 ) and has a kink-like x 3 -dependence, localized at the monopoles. For x 3 — > ±00 
the matrix U is in the Cartan subgroup, see [16] for an explicit example. The asymptotic 
form of the BPS equations (2.10) determines the residual fields to be of the form, 

As^iUdaU- 1 , Ae-iUdeU" 1 , <p « U^U' 1 . (2.12) 

In addition to the BPS equations (2.10) one finds that asymptotically T>z^% ~ T^w^i ~ 
D w 0. 

Discs T)±: On the discs at infinity the fields approach pure vortex behavior exponentially 
fast, with suppressed corrections ~ e -\AmAx 3 \ ^g^g ^fie- distance to the monopoles). 

Therefore one has B\ w B2 ~ -4,3 w A r ~ and tp w 4>q. The nontrivial fields Bs,Aq in 
general take different values at the two discs at infinity, but not the abelian U(l) part. In 
particular one has, 

Ei m b± ~diag(a± ...,a±) , A 9 \v ± = A^ , Tr A e \® + = Tr A e \v_ , (2.13) 

and analogously for £>3|d ± . The BPS equations (2.10) imply then [^4 a =i,2, 4>o\ ~ so that 
[B3, 4>o] ~ 0. In addition one finds P^S, w 0. 

2.3 Fluctuation operators 

The main focus of the investigation in [8] was on the quantum properties of confined 
monopoles. For this purpose the full N = 2 SQCD Lagrangian had to be expanded 
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to second order in the quantum fields in a background that satisfies the BPS equations 
(2.10), i.e. 



Sa = Sji + sn , 



V2<p = V2 (p — 0,4 + ids 

Si2 = Si2 ■ 



(2.14) 



With a convenient gauge fixing term the Lagrangian quadratic in the quantum fields can 
be written in the form 



C 



(2) 
bos 



Tr{W^ (dl + LL r ) W} + ... , 



where the four-component fluctuation field W is given by, 



W 



1 lT 

g a z , Si 1 , Si 2 I 



(2.15) 



(2.16) 



with the complex coordinates z,w defined as in (2.10). We emphasize that here also the 
trivial fluctuations sii>jvy and Si2 around the vanishing backgrounds ^i>N f = £12 = are 
included. This will be discussed below. 

The ellipsis stands for terms of a similar structure. These are ghost fields and the 
trivial 05 and ao fluctuations, which are governed by the one-one matrix component of 
L^L. There are no zero modes from these fields, as we will see. The full supersymmetric 
Lagrangian also includes fermionic fluctuations governed by operators L and L\ which 
contribute zero modes for the classical e.o.m., and are in fact the reason to organize the 
bosonic fluctuations according to (2.15). We refer to [8] for more details, which are not 
relevant for the present considerations. 

The objects of main interest are the fluctuation operators L, which depend on the 
classical background fields defined by the BPS equations (2.10). Their detailed form is 



L 



Lt 



-if 



-ig^l 

TTSf 



where we introduced Euclidean quaternions and the abbreviations, 



and a k = (cj fc ,-il 2 ), t> = a^V k . 



(2.17) 



(2.18) 



The index k runs now over four values and P4 was introduced below (2.10). The super- 
scripts a, f indicate the adjoint and fundamental representation, whereas the superscript r 
indicates action from the right. Explicitly, the right action on an adjoint field X is just ma- 
trix multiplication, ■ X := XEj, and on a fundamental field it is tensor multiplication 
X/ • yi := yi <g) Sj. As usual, summation over repeated flavor indices is implied. 

The operators act in the space of the direct sum of adjoint and fundamental fields, for 

V" V"/ 

example L : y — > y , . It is with regard to the natural scalar product in this space that 

the adjoint operator is the hermitian conjugate of L, and vice versa. 

We want to emphasize that the fluctuation operators (2.17) originate from the quadratic 
Lagrangian in the BPS background and thus describe the fluctuations w.r.t. the full, second 
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order field equations. This is in contrast to the fluctuations of the g BPS equations (2.10). 
We will discuss this subtlety in more detail below. 

In order to compute quantum corrections the quantity of prime interest is the dif- 
ference in the spectral densities of the operators LL^ and L^L. We will call this quantity 
henceforth "spectral density" Ap. This spectral density can be conveniently extracted from 
an index theorem. The advantage of the index theorem calculation is that usually only the 
asymptotic behavior, i.e. the topological properties of the classical background fields, have 
to be known. In the case of confined monopoles this is not a straight forward issue as we 
will discuss now. 

3. Index Theorem 

The index of an operator L , Ind(L) = n^ ti — tl° ll ^ , counts the difference in the number 
of zero modes of 3 L and L\ and can be obtained from an IR-regulated expression: 

IiM > := MlJZTw-TiJTw} =* M(L)=Kmi(M), (3-D 

where the modified symbol for the trace indicates that the trace is taken now also over the 
functional Hilbert space. Of particular interest for us is the application to non-compact 
spaces, as it was developed in [11, 12, 13]. 

The spectral density Ap can be extracted by a Laplace transformation from the con- 
tinuum contribution to the index function X{M): 

/— M 2 
| u}2 + M2 Ap, (3.2) 

where the signs are chosen such that they match the conventions of [8]. The measure for 
the integration J over the continuum-mode eigen-values u is defined via the l.h.s. 

The usual technique to compute the index is to transform I(M 2 ) into an anomaly 
term, which can be evaluated for M —> oo, and a surface term, which can be conveniently 
computed if the fields are asymptotically trivial. The index is therefore determined by 
the topological properties of the background fields. As discussed above, the situation for 
confined monopoles is rather different. In the following we will introduce the necessary 
generalizations of index theorem calculations and develop techniques for the case of con- 
fined monopoles considered here. 

The BPS equations (2.10) imply a certain structure for the fluctuation operators and thus 
the index function X{M) given in (3.1). The basic input for the index calculation are not 

3 A simple argument for the norm of zero modes also shows that iJ Lipo = L%/jo = and LL^tpo = 
O Ltyo = 0. 
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the operators L, itself but their products, which are of the form 



Li) 



ig {TP^iY 
-ig (t^iY 



-ig (#>£;) r 

-Sijpp + g 2 EjZi 

ig (PZjY 



(3.3) 



x' 
y' 



They also have a nontrivial matrix structure in flavor space and act by 

The j-BPS equations (2.10) imply the following relations for the building blocks of these 
operators: 

~-V 2 k -a 3 D , pp = V% + <T k (2B k + 8 kt3 D) , 





(3.4) 



With Tp being the hermitian conjugate matrix, one obtains the remaining combinations for 
the expressions in the second line of (3.4). 

Contrary to the usual situation for |-BPS backgrounds neither of the fluctuation op- 
erators L, (2.17) is necessarily positive definite. This can be seen from the (positive) 
norm of these operators: 

2 



L 


~x~ 











V 



">x\\ + \\t> 



t>x + \\t> 



Vi 



| 2 + g 2 { \\yi ® 

-ig(y im p^X r 



| 2 + 9 2 ( hi ® ^ 



m Ui 



+ II^^H 2 ) 

f^x m n ), 



Vi 



(3.5) 



where the second line in both equations give the cross terms of the norm square which are 
not necessarily positive (we indicated the gauge group indices to, n explicitly). Usually for 
one of the two operators L, this cross term vanishes. In that case the positive definite 
non-derivative terms imply that the l.h.s. vanishes only if the state [X, itself is zero. 
Therefore the respective operator has no nontrivial zero mode. In the case at hand one 
can see from (3.4) that neither for L nor for the non-positive cross terms vanish and 
therefore both operators in general will have non-trivial zero modes. 

3.1 Principle structure. Anomaly vs. anomaly 

Before actually calculating the index function X{M) given in (2.17), we discuss the well 
known basic structure behind such calculations [11, 13, 12]. This will serve as refer- 
ence point for the necessary generalizations for the computation of the index for confined 
monopoles. 

4 We are interested in possible zero modes here. Therefore we can safely neglect surface terms which do 
not contribute for such normalizable and localized discrete states. Non-normalizable zero modes that lead 
to many subtleties for non-abelian monopoles [18] do not alter the following conclusions. 
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First, one introduces the auxiliary "Hamiltonian" which factorizes the original opera- 



tors, 



H 



L 



--: T% + K{x) 



-H 2 + M 2 



L^L + M 2 



Ltf + M 2 



(3.6) 



where we note that (H + M)(—H + M) = —H 2 + M 2 . The matrices Y l are assumed to 
satisfy a Clifford algebra, i.e. {r l ,r J } = 5 13 . Below we will discuss also situations where 
this is not case. The size and number of these matrices is kept unspecified here. The 
relation on the right in (3.6) and the factorization given beneath, allows one to write the 
index function (3.1), which is defined in terms of second order differential operators, as a 
functional of a first order operator. Restricting the trace in (3.1) to the component and 
gauge indices one finds for the local index function J'(x,y, M): 



G M (x,y) : 



1 



H + M 

J(x,y,M) 



(-H + M) 
Tr{r 5 



1 



-H 2 + M 2 

M 2 



H 2 + M 2 



} = Tr{T 5 MG M (x,y)} . (3.7) 



Here we have introduced the Green's function Gm(x,u) and the chirality matrix T§ = 
diag(l, — 1). The symbol Tr includes now also the trace over the gamma matrices. The 
global index function (3.1) is then given by I(M) = J x J{x,y,M)\ y=x , but it will be 
beneficial to keep the local expression. 

Second, one uses the defining equations for the Green's function to rewrite the local 
index function J{x,y,M) as a total derivative. The Green's function Gm(x,v) satisfies 
the following first order equations, 



[H x + M]G M (x,y) 
G M (x,y)[Hl + M] 



[T%i + K x + M] G M (x,y) = 5{x - y)l 
G M (x,y)[-T l d yl +K y + M] =S(x-y)t 



(3.8) 



where we have indicated the coordinate dependence of the operators by an index and 1 on 
the r.h.s. stands for the identity in spinor- and gauge-matrix space. The off-diagonal block 
structure of the operators imply that they anti-commute with the matrix, i.e. 



{r 5 , H} = {T 5 ,T i } = {T 5 ,K} = 0. 



(3.9) 



Therefore, taking the trace of the two equations in (3.8) with T$ and adding the result 
gives, 



J(x, y, M) = -\ (d xi + d yi ) Tr { T 5 G M (x, y) } 

-^Tr{r 5 (K x -K y )G M (x,y)} 



(3.10) 



where we have used the cyclicity of the trace in the finite-dimensional vector spaces, and 
that Tr{ T§ 1 5(x — y) } = 0. This relation is of course problematic in the limit y — > x, as is 
the second line in (3.10) which seems to vanishing in this limit. These ambiguities are the 
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source of possible anomalies and regularization is necessary for a proper treatment. Putting 
this subtlety aside for a moment one sees that the index function Z(M) = f J{x, y, M)\ y=x 
is given by the integral of a total divergence and thus determined by a surface term, i.e. 
by the topological properties of the fields in the operators L and L*. 

It has turned out that Pauli-Villars regularization is most convenient in the given 
context. This amounts to replacing the Green's function given in (3.7) by 

G M (x,y) — > G^(x,y) = G M (x,y)-G fl (x,y), (3.11) 

where the regulator mass fi is sent to infinity at the end. For the Pauli-Villars Green's 
function G^(x,y) the same relation as in (3.10) holds and the difference of these two 
equations yields 

J(M) = / J(x,y,M)\ y=x 

J X 

= lim / J(x,y,fi)\y= x - I lim / diji(x), (3.12) 

^->oc J x ^oo J x 

where the problematic term Tr{ T$ 1 5{x — y) } has canceled. The regularized version of the 
last term in (3.10), which is ~ (K x — K y ) G reg (x, y), rigorously vanishes for y — )■ x since 
G reg (x,x) is well defined for every finite \x. The potential anomaly has been shifted into 
the first term of (3.12). The current ji{x) is given by 

ji(x) = Tr { T 5 I\ G rc §(x, y) } \ y=x . (3.13) 

Finally we want to emphasize a point which has to be generalized in what follows. The 
relation (3.10), which eventually allows one to define the index function as a surface term, 
is obtained in this form thanks to the fact that the masses M and \x are scalars that com- 
mute with all other quantities. 



Anomaly vs. Anomaly. We briefly outline the well known relation between index theo- 
rems and chiral anomalies. The purpose is to emphasize a particular observation concerning 
mass corrections and central charge anomalies for solitons (kinks), vortices and monopoles 
[19, 20, 8]: As has been mentioned, the nontrivial mass corrections and central charge 
anomalies stem from the spectral density (3.2), which is non- vanishing only if the index 
function X{M) is not independent of M. On the other hand, equation (3.12) shows that 
the anomaly part, the first term on the r.h.s., is independent of M. We will argue now 
why the non/existence of a chiral anomaly in the auxiliary model associated with the index 
under consideration, in general implies an anomaly/vanishing corrections for the solitonic 
objects in the actual model. The Euclidean auxiliary model is: 

1 S| ux = y d D x E iTr{&{H + M)*} G M {x,y) = (¥(*) &{y)) , (3.14) 

where GM(x,y) is the Green's function defined in in (3.7), here given in terms of the 
Euclidean two-point function of the auxiliary system. The associated Minkowski action 
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has a chiral symmetry, broken by the regulator mass M, whose current is 5 

jfa = iTr{^TT 5 ^} d^jfa =2^MTV{^r 5 *} -a (5) (x). (3.15) 

The second relation is the anomalous conservation equation, including the explicit breaking 
term, for the chiral current. Taking the expectation value of the anomalous conservation 
equation one sees, after Wick rotation x° = — ix 1 ,^ = — iF , that the first term on the 
r.h.s. is equal to — 2Tr {T§ M Gm} = -2 J(x,i, M), see (3.7). Hence the index function 

(3.12) can alternatively be written as 

X(M)= f J(x,x,M) = -l f a (5) (aO-i / $ (i{ 5) ) ■ (3-16) 

J X J X J X 

Given the definition of the chiral current (3.15) and the propagator in (3.14) one sees that 
the expectation value of the chiral current is identical with the above encountered current 

(3.13) , i.e. after Wick rotation is {jfy) = j l ( x )- As it is well known, the chiral anomaly 
is associated with the index of the Euclidean Dirac operator ij)^ := H of the Lagrangian 

(3.14) . Explicitly one has 

index^) =Z(0) = -i J a (5) (x), (3.17) 

see for example [22] for a recent detailed account. This means that if the auxiliary system 
(3.14) has a chiral anomaly the index function Z(M) is completely determined by it and 
there are no contributions from the surface term 6 in (3.12). As mentioned above, the 
anomaly term is independent of the IR-regulator mass M and consequently the continuum 
contribution (3.2) and thus the spectral density Ap vanishes in this case. 

The quantum mass corrections and central charge anomalies for solitons are determined 
by the spectral density Ap. Hence, the arguments given here confirm the observed fact 
[19, 20, 8], that if the auxiliary system in the soliton background is anomalous the quantum 
corrections for the soliton mass and central charge anomaly vanish, and vice versa. 

In particular, since there are no chiral anomalies in odd dimensions, one finds in general 
nontrivial mass corrections and anomalous central charges for solitons that occupy one and 
three spatial dimensions like kinks and monopoles. An exception to the rule occurs when 
the field content leads to a vanishing overall factor, as it is the case for conformal models 
like M = 4 supersymmetric YM-theory [20] or M = 2 SQCD with 2N = N f [8]. 

3.2 Overall anomaly 

Since the considered setting is three-dimensional the anomaly part in (3.12) vanishes. We 
sketch here the proof of this statement and give a more detailed account in appendix A. 

5 In Euclidean space ^ and ^ are independent, whereas going to Minkowski space implies — > fy, 
see e.g. [21]. The Euclidean and Minkowski space Green's function are related as Gm(x , y , . . .) = 
i GS ink (^> V°, ■ ■ 01 fCZ T y i ■ Hereby is Gf/ nk = (* *} and * = vf^r is the Dirac conjugate spinor. 

6 The only way around this would be to assume that the surface term contribution vanishes only in the 
limit M — > 0. This however leads to a contradiction with (3.2). 
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The anomaly 



A tot = f j,no^ x)= Um f Tr{r 5 jf i A 



2 J ly =a: ' 



(3.18) 



is easily evaluated by expanding the Green's function (-H 2 + /i 2 )" 1 for large ^. The 
notation A tot indicates that this is the total anomaly of the system, below we will find a 
different and essential sub-anomaly. For the given fluctuation operators (2.17) H is of the 
form (3.6) with 



Is 



a 



a 



, K(x) 



Q 



o- ^ 



(3.19) 



The expansion for large \x reads as 



_ n l, 2 = A M (x - y) + I A v( x ~ Zi)0(z 1 )A IM (z 1 -z 2 )... 0{z n )^{z n - y), (3.20) 

with being the three-dimensional propagator (A.l). The first term obviously vanishes 
under trace with and the insertion O is the deviation from the free Laplacian, 



O = H — d = T l diK + K + 2KSi . 



(3.21) 



Here Ki = -{Tj, K}, which is of form diag(a, a) such that Tr {F^K^ . . . Ki n } = for any 
n. In order to estimate which orders in the expansion (3.20) can contribute in the limit 
fj, — > oo we note that the n'th term scales as /i( 1-2n+r ) for y = x, with r being the number 
of derivatives of O that appear in the product 7 . The term with the maximum number of 
derivatives, r = n, vanishes after taking the trace because of the just mentioned property of 
Ki. Hence, the terms that can contribute are (n = 1, r = 0) and (n = 2, r = 1). Both terms 
are regular for y = x and therefore the second term with the single derivative vanishes by 
symmetric integration. Thus the only a priori surviving (and potentially divergent) term 
is 

A* 



J &no (x,x) 



lim 

fj,— >oo 



8tt 



Tr {T 5 rdiK(x)} 



(3.22) 



which vanishes due to the fact that Tr { r 5 P K(x) } = 2 tr ( A\ - A{ ) tr cr* = 0. 

Thus the overall anomaly contribution to the index function I(M) vanishes, A tot = 0, 
as expected from the discussion following (3.14). Consequently, if non- vanishing, I(M) 
will be M-dependent and lead to a non-vanishing spectral density. 



7 Changing the integration variables as Zi — > x + Zi//j, one sees that the coefficient functions of 0(z\) can 
be evaluated at x. An expansion introduces negative powers of ^ while the ^-integrations are finite due to 
the exponential decay of the propagator. 
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3.3 Index of a boundary term 

Having shown that the anomalous contribution vanishes the index function (3.12) reduces 
to the surface term: 



1(M) 



I hm 



jU-XX) J x 



/ diji(x) 

J X 



\ lim 



d6 dx r j r 



--: X 3oo (M)+[X s+ (M)-T :D _(M) 



/ 



+ / d6 rdr js 



(3.23) 



where r,6,x 3 are the cylindrical coordinates. In the last line we introduced a notation in 
accordance with the definitions of the boundary at infinity in figure 1. To evaluate the last 
two terms in (3.23) one has to know the field configuration on the whole discs at infinity 
D±, which describe the full vortex dependence and is not even known analytically We will 
now show how to resolve this situation and how to reduce the computation of the index 
to the topological properties of the field configuration, as it is the case for more familiar 
situations. 

Cylinder Wall 3oo: 

For the given situation this contribution is easy to evaluate. With (x a ) = r (cos 6 , sin 6) 
one has, 



fl— >oo 



lim / dddx 3 x a j 



with j a = Tr {r 5 r Q [ 



1 

H + M 



( M )]}, (3.24) 



where we indicated the subtraction of the same term with M replaced by fi. To compute 
this quantity we first note that with the asymptotic behavior as discussed around the 
equations (2.11) and (2.12) the building blocks (3.4) simplify considerably. Hence the 
operators (3.3) take the simple ultra-diagonal form, 



\ 2 ®[-Vf+g 2 (T li ®t li y] 



(3.25) 



where ~ again stands for equal up to exponentially suppressed terms ~ e gvr . The relevant 
term for computing j a \ r =oo therefore becomes, 



Tr { r 5 r c 



i 



} 



Tr{r 5 r Q (-F)- 



i 



} 



trcr c 







(3.26) 



Consequently there is no contribution from the cylinder wall, i.e. I? loa (M) = 0. 
The Disc contribution at 2)±: 

With the anomaly and the contribution at the cylinder wall 3oo vanishing the index function 
(3.23) is entirely given by the contribution from the infinite discs T>± at x 3 — > ±oo: 



T{M)=Iv + {M)-Z 3 _{M) 



(3.27) 



- 14 - 



where the explicit form of the disc contributions is with (3.23), (3.13): 



2b ± (M) 



.73 



lim 

— ^OO J 



Tr{r 5 r 3 [ 



.73, 



1 



H + M 



(M)]}l 



^• rcg (x,y)| 



J/=X • 



(3.28) 



With the asymptotic behavior at the discs D± as described around (2.13) and the 
properties for the fundamental building blocks (3.4) the operators L^L and LL^ have the 
following block structure, 



* * 

* * 



(3.29) 



Here and in the following stands for equal up to terms that are exponentially suppressed 
as e ~\ AxAm \, see the discussion around (2.13). We therefore perform a transformation such 
that, 



t\2 



-{OHO 1 ) 



* * 



* * 

* * 



(3.30) 



where O is orthogonal 8 . The resulting harmonized block structure is described in terms of 
the following substructure of operators: 



C 



-V% +igfll 
-igY.\ V\ 



-igT,\ 



-vi 



(3.31) 



where from now on it is always understood that all objects are considered to be at the discs 
D±. Another structure which will be important in the following is given by the operator 



M 



'2 r 



V2{4>o + mi) 



(3.32) 



where 4>q = <po T is the vacuum (2.5) which is a matrix in the fundamental representation 
and (4>o + rrii) is the operator when acting on a field yi in the fundamental representation 
which carries a flavor index i. The operator </>q is the associated adjoint representation, i.e. 
it acts as commutator or written as a matrix it is of the form (<Pq) A b = — 



7 

The operators introduced in (3.31) and (3.32) satisfy the following relations: 

[£,M]Ks[£,d 3 ]KsO, 



C fCAB 



(3.33) 



with C) satisfying the same relations and we note that C) £ is diagonal. These relations 
follow from the asymptotic BPS equations given below (2.13) and the fact that the x 3 
dependence is exponentially located at the position of the monopoles, which are assumed 
to be infinitely far away from the boundary. 



The explicit form of the transformation is O 



[I a b ] witha=(Jg8S)and6 = (8S;8) 
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Transforming also the gamma matrices in (3.28) with the orthogonal transformation 
O the current on the discs D± takes the form, 



j(M, x, y) = Tr {r' 5 r' 3 _L_ } = Tr {r' 5 r' 3 (-h> + m) _ H J + M2 } 

= 2Tr{M [ Cic _ d 2 + M 2 + M 2 ~ Ctf-dl+M 2 + Np\ ^ {3M) 

The actual quantity needed to compute Zj> ± (M) is j Teg (x, y) = js(M, x, y) — j^{n, x, y), see 
(3.28). 

The quantity (3.34) reminds also of our starting point, the definition of the index 
function (3.1). It is this similarity that we are going to use to compute (3.34). The 
difference to (3.1) is that in addition to the parameter M it contains the operator M 
and especially the derivative operator 83, which contrary to the operator C, does not act 
in the bulk of the discs 2)±, but orthogonal to them. Besides these subtleties, which we 
address in a moment, the quantities j(M, x, y) may be considered as the definition of a new 
index function, analogous to (3.1) or more accurately analogous to the local index function 
J{x, y, M) in (3.7), and we will compute it in this spirit. 

We first recall the basic properties that were used in section 3.1 to formulate the index 
function of the form (3.1) as the sum of an anomaly and surface term: i.) The factorization 
property which allows to express the index function in terms of a Green's function of a 
first order operator (3.7) which satisfies equations of the form (3.8), and ii.) commuting 
properties with a chiral matrix T§ (3.9) which lead to (3.10). In this last step it is also 
important that M commutes with all other expressions. 

In a first step one has to deal with the presence of the unusual derivative operator 
$3 in (3.34) which acts orthogonal to the boundary and prevents a proper factorization 
mentioned under i). To this end we introduce the Fourier transformed current j in the 
following way: 

j(M,x,y)= J l^itK,//) with 

](M„) : = 2Tr{M[ £t£ + M 2 + M£ ~ C0 + M* + M* ] > ' ^ 

where in the first line we indicated explicitly the coordinate dependence with a = 1 , 2 
denoting the coordinates on the discs D±. From now on we will suppress this notation, as 
in the second line of (3.35). The momentum dependent IR regulator is Mp = M 2 +p 2 . 
Using a complete eigen system of the (diagonal) operator C) C and the fact that it is 
isospectral to the operator ££< (zero-modes are treated without problems separately) it is 
easy to prove the identity (3.35). 

We follow now the steps of section 3.1 which allowed one to express the original local 
index function J{x, y, M) as a surface term and an anomaly. First one has to find the 
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proper factorization in terms of an auxiliary Hamiltonian. Defining the operators 



H 



H 



1 2 ® M + o- 3 ® iM p l 2 -la 8) £ f 

1 2 ® £ 1 2 ® - o- 3 ® iM p l 2 



-1 2 <g>.M + cr 3 <g>iM p l 2 -l 2 <g>£t 

1 2 ®£ -1 2 <8> M - a 3 iM„t 2 



(3.36) 



one can rewrite the Fourier transformed current (3.35) in the desired form: 

j( Mp )=Tv{T 5 n-t—}=-Tr{T 5 G Mp } with: G Mp := I 
— ri ri H. 



(3.37) 



This is analogous to (3.7). We now have to find the generalization of the steps which 
led from (3.8) to (3.10) for j(M p ) and j(n p ) such that we eventually obtain a regularized 
expression analogous to (3.12). The auxiliary Hamiltonian is of the form 

H = T a d Q + K(x) + K Mp with K Mp = -1 4 <g> M + a 3 ® cj 3 ® iM p l 2 , (3.38) 

whereas 9 T", K{x) are independent of and M p and are block off-diagonal so that they 
anti-commute with T^. The operator H is of the same form with Km p — > ~Km p - 

Second, we again use the defining equations for the Green's function to rewrite now 
j(M p ) as a surface term. The Green's function (3.37) satisfy the equations 



[F a d xa +K X + K Mp } G Mp (x, y) = S(x - y)l 
G Mp (x, y) l-r a 8 ya +K y + K Mp } = 8(x- y)t 



(3.39) 



which is similar to (3.8). However, the simple IR regulator M has been replaced by operator 
Km p - Above we took the trace with T$ of (3.8) and the sum led to (3.10). Due to the 
presence of the operator Km p one has to generalize this procedure as follows. We assume 
the existence of a block diagonal matrix T (given explicitly below) such that it commutes 
with Km v and has similar commutation properties as before (3.9), i.e. 10 



{T,T a } = {T,K} = 0=[T,K Mp } 



(3.40) 



The vanishing of the last commutator is the analog of [T5 , M] = in the previous deriva- 
tion. Taking the trace with T of the sum of the two equations in (3.39) one obtains 



Tr {TK Mp G Mp } = ~ U d *°> + d v a ) Tr { rr « G Mp ] 

- \Ti { T{K X - K y ) G Mp } + 6(x - y) Tr T 



(3.41) 



These matrices, though of the same block structure, differ from the ones of section 3.2. They are defined 
via (3.36). 

10 Actually, this relations are stronger than really needed. It suffice that they hold when inserted in the 
trace of Gm„ 
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The l.h.s. should now give j(M p ) from (3.35). To this end we set 



r = 



<g> M 



1 



TT{TK Mp G Mp }=j(M p ), 



(3.42) 



M 2 + M| 



which can be easily seen using (3.37). Because of relations (3.33) this T satisfies the assumed 
commutation relations (3.40). An important property of the nontrivial factor in T is that 
it exists even in the limit M p — > 0, since the zero-eigen values of Ai 2 are projected out by 
the nominator Ai. 

The quantity which is actually needed is not j(M p ) but the regularized expression 
jreg _ j(M p ) — j(fJ-p), see (3.28). One finds the same relation (3.41) for Gp, by replacing 
M p by fi p in the above derivation. The last relation in (3.42) does however not hold in this 
case due to the form of T. Subtracting from (3.41) the analogous equation for Gp p one 
finds 



Tt{T(K Mp Gm p -K^G^)} = - l(d xa +d ya )Tr{TT a (G Mp - G, p )} 

- ^Tr { T (K x - K y )(G Mp - G, p )} . (3.43) 



The terms on the r.h.s. are regular for y — > x since only the regularized Green's function 
G reg = Gm p — Gf_i p appears. For this to be the case it is important to use the same matrix V 
for both cases of the relation (3.41), either with parameter M p or fi p . The second term on 
the r.h.s. of (3.43) vanishes therefore in the limit y — > x and will be not considered further. 

The l.h.s. of (3.43) is not yet of the desired form to give j reg and also it seems that the 
usual anomaly term is missing. Using the explicit form (3.38) of Kx=m p ,^ p and (3.37) the 
l.h.s. of (3.43) can be written as: 



Tr { r (K Mp G Mp - K^Gp,)} = j reg - TV { 1 4 ® ^ T 5 G, p } . (3.44) 



Consequently, one can rewrite the surface term contribution to the index function 
I(M), which comes exclusively from the discs l D±, as a surface term at the boundary of 
the discs dD±, where all fields assume values in the vacuum moduli space plus an anomaly 
term: 



q a (x) = -iTr {TT a G^}\ y=x , a(x) = Tr{t 4 ® *~ T 5 G flp }\ y=x . (3.45) 



We state here only the results for the evaluation of the expressions in (3.45). The proofs 
for the following statements are given in appendix B. Both objects, q a {x) and &(x), have a 
prefactor A4(A4 2 + M^)~ l inside the trace 11 which is of utmost importance and represents 

11 For q a (x) it is explicitly contained in F (3.42). For a(x) it follows from the first equation in (3.35) in 
combination with the already present prefactor, see (B.5). 




with: 
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the modification of the expressions (3.7) and (3.13) that enter in standard formula for the 
index function I(M) (3.12). This factor renders the p-integration finite and projects out 
any poles and massless fields in the propagators that might appear in the M — > limit, 
see the comments below (3.42) and appendix B. 

It turns out that only the anomaly contributes in (3.45) and the result is 

1(M) = A disc = ~ [ Tr( , = = B 3 x r* }, (3.46) 

where we used that X = — see (3.27). Remarkably, and for the non- vanishing 
of the spectral density (3.2) essential, this anomaly is not independent of the IR-regulator 
M. 

3.4 Topological charges 

We describe now the relation to topological quantities for the index (3.46). 
Topological charges: 

In order to define define certain topological invariants we consider the classical energy of the 
field configurations considered here. The BPS equations (2.10) imply that the (classical) 
energy of confined monopoles is given in terms of the total vortex number and the magnetic 
charge of the monopole [8] : 

E cl = 2irv 2 k mv L + j, J d 3 x 8 k Tr {^2 <p B k } , (3.47) 

where the first term is the total vortex tension times the regulated extent of size L in the 
x 3 direction. The total vortex number is k vor = ^- J d 2 xTrBs = Trw, see (2.12). However, 
in the following the second term, which gives the energy of the monopole and encodes the 
magnetic charges will be of prime interest. 

For the given asymptotic behavior (2.12) the surface term in (3.47) reduces to con- 
tributions from magnetic flux through the discs D±. Following [14] one can express this 
flux in terms of the individual contributions to the vortex number according to the sym- 
metry breaking pattern (2.7). Corresponding to the group Hc+f (2-8) there are q distinct 
topological quantum numbers 

1 f q 
A^ or = _ / Tr {£ 3 t° r } V Ak; or = k mr - k mr = , (3.48) 

where t® = diag(0, . . . , l nr , ... 0) is the U(l) generator of the r'th factor of Hc+f (2.8). 
The vacuum value of the adjoint scalar (2.5) can thus be written as <fo = —Ylr=i rnr ^r- 
The monopole mass contribution to the classical energy is then given by 

M c ? on = Jr [ Tr {^2 O S 3 } = -^V^m.Afcr, (3.49) 



The above discussion relates the (multi) monopole mass to a set of topological charges 
which are specific for the case of confined monopoles. There is a more generic set of 
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topological charges which can be assigned to Coulomb (unconfined) and Higgs (confined) 
monopoles. The mass is determined by the behavior of the fields at the boundary at infinity 
which we denote by doo. For the monopoles in the two different phases the boundary is of 
the form 

_(Sl ...Coulomb 
d °°-\D + US)_U3oo .-.Higgs ' (3 - 5 ° j 

where the first line denotes the sphere at infinity for the asymptotic Coulomb monopoles 
and for the confined monopole the boundary 3oo gives no contributions (3.49). At the 
boundary <9oo the fields commute, i.e. [Bk , ^Soc = 0, see (2.13) and [23, 24] for the 
Coulomb monopole. The Cartan subalgebra can be therefore chosen such that it contains 
both asymptotic fields. 

The scalar field at the boundary is V^ipg^ = y/2(f)o = h- H + 0(r~ 1 ), where H = (Hi) 
are the Cartan generators and h is a constant vector 12 . The square root y/2 is a convention 
to suite the definition (3.32). The asymptotic chromo-magnetic fields are of the form: 

Coulomb: B 9oc -> B r = — j-^ g ■ H + 0(r" 3 ) , 

Higgs: B doo -> Bf = ^d r Af at 9± . (3.51) 

In the Coulomb case r denotes the spherical radial coordinate and g is a constant vector [25] 
called "magnetic weight" [23] . The form at the boundary for the confined monopole follows 
from (2.13), where in this case r is now the cylindrical radial coordinates. The gauge fields 
lie also in the Cartan subalgebra and we associate to them a magnetic weight via their 
constant values at the boundary of the discs D± (2.12) as follows: Ae\d® + — Ae\d®_ =: 
2^ g • H. We note that the J7(l)-part of the U(N) gauge field drops out of this definition, 
as it does in (3.48). 

The magnetic weight g is associated with the non-abelian magnetic flux through the 
surface at infinity, 

Foo = f da k B k = i f g H = f B 3 = g • H , (3.52) 



9oo 47r J Si, JV+-Z)- 

where we have assumed that the weights coincide for the Coulomb and Higgs monopole. 
Note that in both cases the flux Foo contains no contribution from the U(l) factor (for 
Coulomb monopoles it decouples from the beginning). Thus if working with the U(N) 
Cartan subalgebra, see appendix C, the weights have to satisfy in both cases Yli9i = 0- 

Thus with every confined multimonopole in the Higgs phase one can associate a 
Coulomb multimonopole by identifying their magnetic weights and thus their total flux 
to infinity. They have then equal masses: 

g h = g c =: g => M h mon = M c mon = \h-g, (3.53) 



12 The approximation to order 0(r _1 ), r being the spherical radial coordinate, applies to Coulomb 
monopoles [14] and is the maximum necessary information needed. For the confined monopoles the asymp- 
totic value is reached with exponential precision at the relevant boundary ®± (2.13). 
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where we assumed that for both cases the vacua <f>o are identical, for example (2.5), which 
gives h = (hi = — 2mj). The first term in the energy (3.47), which is proportional to the 
spatial extent of the system and is present only for confined monopoles, carries the infinite 
energy of the confinement transition and is exclusively given by the U(l) factor of the 
gauge group. 

In [23] it was shown for Coulomb monopoles that the magnetic weight has to satisfy a 
quantization condition in terms of the dual simple roots 



g = 4vr ]T n a (3* {a) with n„eZ and 0? a) = = /3 (a) , (3.54) 

0=1 P(a) 



where (3r a \ are the simple roots, we refer to appendix C for our conventions, and the 
integers n a are the so called GNO charges. If the unbroken gauge group is non-abelian, 
as we assume in general, not all simple roots will have a non-vanishing scalar product 
with h. However, we can choose it to be non-negative, as it is the case for the vacuum 
(2.5): h ■ (3 {a) = V2(m a+1 - m a ) > 0. The GNO char ges are separated accordingly into 
topological and non-topological charges [24]: 



N-l 



h g = 4vr 2^ n a h ■ (3 (a) = in^nth- (5 {t) , (3.55) 



a=l t 



where (3u\ are those simple roots that have a non- vanishing scalar product with h. For the 
vacuum <pQ (2.5) with q groups of n r degenerate masses, see (2.8), there are q — 1 topological 
charges and associated simple roots: 



r 



( n t,P(t)) = (tV,/3(a r )) with a r = ^2n s , r = l,...,q-l. (3.56) 

s=l 

Inserting the total flux relation (3.52) into the definition of the principal topological charges 
for the confined monopole (3.48) one finds a relation with the topological GNO charges if 
one assumes that the magnetic weight for the confined monopole is of the form (3.54), or 
equivalently that the masses coincide 13 (3.53): 

r 9-1 <? 

AC = V-V-! & V=E A C n tot :=^n ar = -^rAC r , (3-57) 

s=l r=l r=l 

where the last relation gives the total magnetic charge. This allows one to associate to 
every confined multimonopole a Coulomb multimonopole, and vice versa, in a unique way. 
For example, a single monopole with one GNO charge, n af = 1 for r = r and otherwise 
zero, corresponds to the vortex charges AkJ OT = 1 = — A/c™^. And a confined monopole 
with minimum set of vortex charges AkJ or = 1 = — AA;™- that obey (3.48), correspond to 
s unit-GNO charges ti . = . . .n af g _ 1 = 1 with total magnetic charge n tot = s. Hence, 



13 By the direct comparison of the masses [14] arrived at similar relations, though there seem to be some 
typos in [14] . To keep the formulas compact we set the non-existing topological GNO charges n ao , n ajv 
(3.56) equal to zero. 
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for a confined monopole with unit total magnetic charge the confining vortices have to sit 
in neighboring group factors of the unbroken group (2.8). This fits the picture that the 
confined monopoles correspond to kinks that connect neighboring vacua in CP n models, 
as it was recently confirmed also at the quantum level [8]. 

Index and spectral density: 

The Index function X(M) (3.46) has two competing contributions, one from the adjoint 
sector and one from the fundamental sector. The contribution from the fundamental sector, 
see (3.32), is easily obtained with the explicit matrix realization (2.5): 



l f (M) 



, N f q 

EE 



V2( 



in; 



i=l r- 



\ y/2(mi 



m r 



+ M 2 



Akl 



(3.58) 



where we used the definition (3.48). Note that in the limit M — > the sum for i > N 
reduces to 14 ~ J2r=i A K° r = by (3.48). Hence, there are no contributions to the actual 
index 1(0) from i > N, but there will be for the spectral density in (3.2). 

The contribution to X(M) from the adjoint sector is obtained by using the relations of 
appendix C. In particular, with \/2<Ao = ^ ' (C6), one has for the adjoint corner in 
(3.32), 



M 



VM 2 + M 2 







iA 
-iA 



with [A] 



(h • a) 



^/(h-cx) 2 + M 2 



(3.59) 



where A is the diagonal matrix where each entry (3.59) corresponds to one of the positive 
roots a, see (C.6) and (C.3). The integration in (3.46) for the adjoint field B| gives the 
non-abelian flux (3.52) in the adjoint representation, i.e. g- H^ d , which is of the same form 
as (3.59) where the matrix elements of A are now given by (g-a). The adjoint contribution 
to the index function I(M) (3.46) is then: 



X a (M) 



1 

"2^ 



E 



(h-a)(g-a) 
y/(h ■ a) 2 + M 2 



(3.60) 



where the sum runs over all positive roots a m <n, see below (C.3). Before elaborating 
on this expression we note that it can be written in a similar form as the fundamental 
contribution. One has, 



h-a H = V2 (m k - m £ ) , g • a ke = 2vr [ (n fc - n fe _i) - (n £ - n £ _i) ] 



(3.61) 



where we used the components of h given below (3.53) and (3.54). In the second expression 
appear all GNO charges, also the non-topological ones. Inserting these expressions one finds 



N 



l a (M) 



EE 



V2( 



rn; 



i=l r- 



\ ^/2(mi 



m r 



+ M 2 



Akl 



(3.62) 



4 We recall the ordering mi < . . . < rriN s of the masses, see below (2.1). 
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where we have used that (n s — n s _i) = n Qr — n ctr _ 1 = A/c™ r . This is except for the 

s=a r _i+l 

difference in range of the first summation, which however does affect the index 1(0) itself, 
(—2) times the fundamental contribution. 

The total index function X(M) can therefore be written as 

AM) = -|( 2 E-E)E AC • (3-63) 

This is the form that was used in [8] to extract the spectral density according (3.2) by a 
Laplace transformation: 

A, = ( 2 £ - fj ) ± J * ^::j 2 AC /(«»), (3.64, 

i=l i=l r=l 

with w 2 = fc 2 + 2(mj — m r ) 2 and f(oj 2 ) an arbitrary function. This spectral density 
determines the full perturbative quantum energies of multiple confined monopoles and in 
addition an anomaly in the associated central charge, notably the magnetic charge in (3.47) 
[8]. As can be seen from (3.63), the fundamental and the adjoint sector compete in their 
contribution to the index. However, this relative sign in the two sums of (3.63) is of utmost 
importance in the quantum theory since it produces in the presence of confined monopoles 
the /3-function coefficient for the coupling constant renormalization. We refer to [8] for 
further details. 

We can now turn to the computation of the index and its relation to the above discussed 
topological charges. For this we consider the expression (3.60). As already mentioned, the 
terms with i > N in X f (3.58) do not contribute for M = and thus up to a factor it is 
identical to the adjoint contribution (3.60). Concretely one has, 

X : = 1(0) = X a (0)+X f (0) = ^X a (0), (3.65) 

where X = Ind(X) = n^ L — n ^ LL ^ (3-1)- Thus the fundamental contribution halves the 
index of the adjoint sector. It is worth mentioning that the adjoint index (3.60) for the 
confined monopoles is exactly the same 15 as for the associated Coulomb monopoles [24], 
which have the same magnetic weight g. It is difficult to say at this point what the halving 
of the index for the confined monopoles means. As it was discussed above, neither L nor 
are strictly positive (3.5), and these fluctuation operators describe the fluctuations of 
the full second order field equations and not just the \ BPS equations (2.10). However, for 
so far in the literature considered g BPS equations these two sorts of fluctuation operators 
were identical. We address this question in more detail in the next section. 

In the evaluation of X a (0) one has to account for the fact that the factors (h ■ a) 
are non-negative, but that they cancel only if they are non-zero, otherwise there is no 



3 The overall minus sign is a convention that we inherited here from [8]. 
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contribution from the respective root. Therefore the sum is restricted to those positive 
roots OL m<n for which {h ■ at) 7^ 0: 



1 



q— 1 a r N 



X a (0) 



Yl ' (9 ■ a) = Y Y Y fa ■ <w) 





(a r + a r -i) Aky 



(3.66) 



r=l 



where the prime in the first sum indicates the discussed restriction. The function a r was 
defined in (3.56) and n r are again the number of degenerated masses in the r = 1, . . . q 
groups. The last line gives the index in terms of the two different topological charges that 
were introduced above. For the case of abelian monopoles, i.e. maximal breaking of the 
gauge group (n r = 1 for all r = 1, . . . q = N), (3.66) reduces to the known result [24], 



where n tot is the total magnetic charge (3.57). 

In [14] some examples of the (framed) moduli space of Coulomb and confined monopoles 
were explicitly constructed via the rational map and moduli matrix construction, and it 
was shown that they are identical. These moduli spaces are, 



AV=D = R3 x Sl * (C^"'^ 1 x CP"'^- 1 ) , M ( i,i )s ep = (C x C*) 2 . (3.68) 



The first one describes the moduli space of a Coulomb/Higgs monopole in the generic 
symmetry breaking pattern (2.8) with unit charge n a . = 1. The second moduli space 
refers to a charge two Coulomb monopole in the completely broken phase of SU(3) or, 
respectively, to a charge two confined monopole with symmetry breaking pattern U (3) — > 
U(l) 3 with ni = ri2 = 1. The index sep indicates that the two monopoles are confined 
by independent flux tubes, in contrast two freely moving monopoles aligned on a single 
confining vortex 16 (in figure 1 this corresponds to the situation with one "dot" of charge one 
at each flux tube). A further example that is discussed in [14] is the charge two monopole 
in minimal breaking SU(3) SU{2) x U(l) or U(3) U{2) x U(l) for the Coulomb and 
Higgs case respectively. 

In all these explicit cases one finds for the (real) dimension of the moduli space 



16 In [14] also the situation ol two aligned monopoles is considered. The moduli space is ol dimension 
six in this case. This does not fit in the the counting that we will propose in a moment. Apparently, the 
boundary conditions formulated in section 2.2 do not account for such solutions where the monopoles on a 
given flux tube can be arbitrarily separated from each other, but describe only multimonopoles which are 
close or on top of each other if they are situated in the same confining vortex. 




(3.67) 



dim(ACod) = dim(AC d ) = 2|X a (0)|=4|X 



(3.69) 



with X a (0) given by (3.66). 
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We are thus led to the conjecture, that in general the index X (3.65) that we computed 
here, gives a quarter of the dimension of the moduli space of non-abelian confined multi- 
monopoles, where on each separated flux tube the monopoles are confined to be close or 
on top of each other, see footnote 16. 

3.5 Parameter counting 

We discuss now the relation of the index which was computed in the previous section to 
the dimension of the moduli space, given by the zero mode fluctuations around the BPS 
equations (2.10). 

The operators (2.17) were defined via the expansion of the Lagrangian around solutions 
of the BPS equations (2.10), and it is the supersymmetry structure, i.e. when taking 
into account the fermionic fluctuations, that determined how the fluctuations (2.16) were 
organized and thus the form of the fluctuation operators (2.17). We have to refer to [8] for 
details. One central result of this investigation is the spectral density (3.64), obtained from 
the index function X{M) (3.63). However, for the evaluation of the index X = X(0) itself 
we showed that the trivial fluctuation with i > N, see (2.16), do not contribute. Therefore 
for the following comparison one has to consider only the first N flavors, or equivalently 
set Nf = N. 

BPS fluctuations 

Even for Nf = N the fluctuations (2.16) still contain the trivial fluctuations Si2 which do 
not appear in the fluctuations of the BPS equations. The fluctuation equations for the 
BPS equations are obtained by inserting an expansion as given in (2.14) for the classical 
fields in (2.10) itself. However, the hermiticity condition (2.9) for BPS solutions implies 
05 = (it was irrelevant also before), and of course Si2 is not present and we call the 
relevant fluctuation sn = s, in the following. In addition one has to implement the gauge 
fixing condition on the fluctuations. The gauge used in [8] to obtain the structure (2.15), 
amounts to imposing the condition, 

Vka k + ^ (si ® % - ^ <g> a,) = , (3.70) 

where the notation k was introduced in (2.18). The resulting fluctuation equations are 
given by, 

£(i) = T^zO-z + T>wO, w + ig 2 Si ® E, = , 

£(2) = T> z Si-ia z Y>i = 0, £( 2 ') = V w Si - ia w Y,i = 0, 

£(3) = *D w a z - V z a w = . (3-71) 

Hereby is Em the fluctuation equation for k = 1 in (2.10) and £^ is a complex combination 
of the k = 1 and k = 2 fluctuation equations. As mentioned below (2.10), the k = 1,2 
BPS equations are the integrability conditions for the two last ones. Similarly one finds 
that fluctuation equations £(2) 5 £(2') imply the fluctuation equation £r 3 y. V w £( 2 ) —T^z£{2') = 
i£(3)£». 
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Therefore, only the first three equations in (3.71) are independent and define the 
fluctuation operator: 



W 



\a z , Si) T 



V 



-v% -vi 

igZl 
ig^l 



-vi 



(3.72) 



where the fluctuation field W is the same as in (2.16) with the last entry, the trivial 
fluctuation Sj2 ; deleted. The ordering and notation for the fluctuation operator will become 
clear in a moment. The operator adjoint to the one in (3.72) reads as, 



L 



VI 



-igZ* 
-igt\ 
£> f - V { - 



1,2. 



Lt : 


' r" 


— > 


~ x'~ 


and L : 


' x ' 




' x' a ~ 




. v * _ 















(2.18). 



(3.73) 



with 

These spaces are all of 



We used here the same notation as beneath 
dimension 3N 2 = 2 d a d + N x df un = d a( j + 2 A'" x df nQ , where (i a d, df un denote the dimension 
of the respective representations of U(N). 

Contrary to the previous situation, one of the two operators is now strictly positive, 



L 


' x~ 











\ 2 +9 2 ( 



\yf ®%f + \\xi: i \\ 2 ) =0 



0, 



(3.74) 



and therefore the index ind(L^) computes the number of zero modes of L* , i.e. the (complex) 
dimension of the moduli space of solutions of the BPS equations (2.10). However, there is 
a relation to the previously considered operators (2.17). Using the explicit matrices (2.18), 
one finds that by deleting certain rows and columns from the operators (2.17) that, 



Lt = L t| 



(2,4) 



L = L\ 



(4,2) , 



(3.75) 



where we have indicated which (row, column) has to be deleted to give the BPS fluctuation 
operators (3.72), (3.73). 

For the first operator L\ the deletion described in (3.75), amounts to setting the last 
component in W (2.16), i.e. the trivial fluctuation Sj2, to zero. The second row of L reduces 
then to the integrability condition £r s \ (3.71), which is implied by the other equations for 
zero modes. The question is if this is consistent, i.e. if all zero modes of are of this form. 
Inspection of the product LL\ which has the same zero modes as U , shows that with 
(3.4) the last component of W decouples. The zero mode equation of for this decoupled 
mode implies Si2 <8> £j = so that there are no such nontrivial zero modes. Consequently, 
the zero modes of and L* are identical, i.e. 



L*W = => W = [X 1 , X 2 , y u Of 



°(Lt)=„0(Lt )j 



(3.76) 



where n° denotes the respective number of zero modes. 
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The situation is different for the second set of operators in (3.75). Inspection of the 
product L^L (3.3) shows that in this case the first component of the fluctuation field 
decouples and that the zero mode equation of L implies X = and thus there is no 
such zero mode. Nontrivial zero modes are thus of the form W = [0, X, yj, yf] T - Again, 
one of the component zero mode equations of L is implied by the other three via the 
integrability condition. The residual three equations can be encoded in an operator 17 L: 



LW = => L 



igYS -V% 
igt\ -VI 



(3.77) 



The deletion of row and column as given in (3.75), on the other hand, amounts to 
setting the second component in the fluctuation field to zero, instead of the first one, 
and in addition one needs a nontrivial equation which is not implied by an integrability 
condition. Clearly the zero modes of L and L are different, the latter one does not have any 
(3.74) whereas L is not necessarily positive (3.5), and according to the conjecture (3.69) 
the number of zero modes is, 

n°(L) = n°{L) = ^n°(L t ). (3.78) 

The operators U and L have a similar structure, with the roles of fundamental and 
adjoint representation interchanged, and a complex conjugation of the entries. One could 
try to find a 2 : 1 map between the zero mode solutions of L' and L to proof the last 
equality in (3.78). Certainly one could also try to compute directly the index of D . In this 
regard we want to mention one difference to the considerations in the previous sections. 
The auxiliary Hamiltonian (3.6) is now with (3.72), (3.73) of the form, 

H = Pdi + K with { P, f j } ^ 2 5 ij , (3.79) 

i.e. the P-matrices do not satisfy a Clifford algebra. A number of manipulations of the 
previous section need only that these matrices anti-commute with T§, which is still the 
case, but many other details will be different. We have to leave both of the mentioned 
considerations for a separate investigation. 



4. Summary and Conclusions 

In this paper we computed the index and the associated spectral density, which plays 
a central role in quantum computations [8], for fluctuation operators which are defined 
via the Lagrangian in a confined multimonopole background. These confined monopole 
backgrounds describe asymptotically nontrivial field configurations. To this end it was 
necessary to generalize the standard index calculations of [11, 12, 13] appropriately. The 
general strategy for such computations is as follows: 

17 In order to emphasize the similarities with U we ordered the fluctuation field as W = [— yf, y\, X] T 
for the definition of L. 
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We assume that the boundary of the space occupied by the background is the sum of 
two components, dM. = 9M tr i v U(9M non t r i v , where on one component the field configurations 
are trivial but not on the other. Then i.) reformulate the index as a sum of an anomaly 
and a surface term using the standard techniques, ii.) The anomaly and the surface term 
on the boundary <9M tr i v can be computed in the limits as outlined in [11, 12, 13]. iii.) 
Rewrite the nontrivial surface term on <9M nontr i v as a generalized index, as in (3.34). iv.) 
Fourier transform the derivative operators which transverse to the boundary 5M non triv v -) 
Find a proper factorization of the auxiliary Hamiltonian for the index on the boundary 
and reformulate this index as the sum of an anomaly and a surface term. If the fields are 
trivial at the boundary of cM non triv) which has to be the case for the above splitting of dM 
to be consistent, all terms can be computed in the convenient limits. 

The resulting index for the confined monopoles is exclusively given by the anomaly 
^jdisc (3_4g) on th e di scs which represent the boundary <9M nontr i v in this case. Con- 
trary to standard index calculations, this anomaly depends nontrivially on the IR regulator 
mass M and thus leads to a non- vanishing spectral density (3.64). We were able to express 
this index in terms of topological charges of the confined monopoles, which we also related 
to topological charges of an associated Coulomb monopole (3.66). From some existing ex- 
amples of the moduli space of confined monopoles we conjectured that the index presented 
here counts a quarter of the dimension of the moduli space (3.69). Finally we compared the 
fluctuation operators defined via the Lagrangian with those of the BPS equations. How- 
ever, a detailed proof of the conjecture regarding the parameter counting has to be left for 
a separate investigation. 

Acknowledgments: I thank A. Rebhan and D. Burke for many useful comments on the 
script. This work is supported by the Agence Nationale de la Recherche (ANR). 



A. Overall Anomaly 

Here and for other proofs we will need some properties of the free massive Euclidean Green's 
function, which we henceforth call "propagator". In D > 1 dimensions the Euclidean 
propagator is given by 

dPk e ik ^ x -y) 



i r d u k e lk{ - x - y > 

A^x-y) :- _ Q2 + ^ - J J2^D^T^2 



1 nU+l) ( " Vk v M*-v\), u=^, (A.l) 



y/n (2tt 2 ) v \\x-y 

where K v {z) is the modified Bessel function of 2'nd kind and for D = 2, 3, 4, . . . one has 
v = 0,5,1,... We note that A /i (Az) = \~ 2u A\^(z). The asymptotic behavior of the 
propagator is dominated by the exponential decay of the Bessel function. The large mass 
limit is of the form 

A(x-v) ^ J_ r(^ + ^) /^-v 2 M | x _ y | \ 1 + 0( i\] (A2 ) 
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The folded product of n + 1 propagators, 

A^ +l \x-y):= [ A tl (x-zi)A ll (z 1 -z 2 )...A ll (z n -y) 

J Z-] ,.,.Zn. 



d D k e ik ^-y) 1 T(n-u) 2(u _ n) 

(2tt) d (k 2 + /x 2 )^ 1 (4^ r(n + l)^ ' 1 ' 

is regular at y = x for n > v = 0, i, 1 . . . and scales as ~ /j 2 ^^^ in this case. The insertion 
of r = 1, 2, . . . , n + 1 derivatives . . . di r in the folded product (A. 3), distributing them 
arbitrarily while acting always on the first argument of the respective propagator, 18 gives 

d D k i r k u . . . ku e ik{ - x ~y) 



a « i K lx . . . K lr e _ A (n+i f _ a 



r 



t t f°° dkk°- l + 



Regularity at y = x is given for 2n — r > 2v = 0, 1, 2 . . . and the resulting expression 
scales as ~ n 2u ~( 2n ~ r \ In this case the integration over the unit sphere SI z?— l = ki/k) 
vanishes for odd r by symmetric integration. 

In computing the overall anomaly in (3.18) one has to evaluate the different contributions 
in the expansion (3.20) for large [i up to orders 0(/U -3 ). To this end one changes the 
integration 19 variables to Zi — > x + Zi/fi. The expansion (3.20) then reads as 



M 2(„-n) 



]T f Axi-z^Oix + f)A 1 (z 1 -z 2 )... 0(x + ^A^zn + /x(x - y)) , (A.5) 



where we used the scaling described below (A.l) and v = ^ here. The operator insertions 
(3.21) are of the form 

0(x + f ) = C(x + &) + fl K h (x + f) d z u , (A.6) 

so that the n'th term in the expansion (A.5) for y = x is of order ~ ^ 2 ^^( 2n ~ r ) with r 
being the number of derivatives of (A.6) that appear in the product. The expansion of 
the coefficient functions in (A.6) around x gives further inverse powers of fi, while the 
z-integrations for these terms are finite due to the exponential decay of the propagators. 

As mentioned in the main text, for any given order n in the expansion (A.5) the term 
with exclusively derivative insertions vanishes under the trace with IV Thus in three 
dimensions {2u = 1) only the terms with (n = l,r = 0) and (n = 2,r = 1) survive the 
/x — > oo limit. Both terms are regular 20 at y = x. 



18 Due to translation invariance one may always assume that the derivatives act on the first argument of 
the propagators in (A. 3). 

19 We follow here [11] with the difference that the insertions O contain also derivative operators. 

20 Only the (n — 1, r = 1) term diverges at y — x in the expansion (A.5), but as noted its trace with Ts 
vanishes. 
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Noting that Tr{r 5 AT 2 } = 0, see (3.21) and (3.9), the (n = l,r = 0) contribution to 
the anomaly can be written as 



Jff.°(x,x) = lim (jl [Tr{T 5 T i d i K(x + l)}A 1 (-z)Ai(z) 



lim 

[i— >oo 



JL 

8vr 



Tr {T 5 rdiK(x)} 



(A.7) 



In the second step we expanded the insertion around x and used (A. 3) for v = \. The 
order ~ zj ' [i term in the expansion of the insertion would survive the \x — > oo limit but it 
vanishes by symmetry of the ^-integration. 

For the (re = 2, r = 1) term the coefficient functions of (A. 6) can be directly evaluated 
at x. The resulting term is proportional to a folded product of three propagators with 
a single derivative insertion and thus vanishes by symmetric integration, see (A. 4). This 
leaves for the overall anomaly the contribution (A.7), which vanishes due to the trace with 
r 5 , see (3.22). 

B. Disc anomaly and Surface term 

As shown in (3.45), the index function I(M) for confined monopoles reduces to a surface 
term and anomaly contribution on the discs at infinity %)±. We give here some details in 
the evaluation of these contributions. 

In what follows one needs the explicit form of the products CC\ C^C on the discs T)± 
(3.31). Both products contain an adjoint and fundamental "mass term" which can be seen 
in the diagonal elements of the original operators (3.3). With the behavior at the discs D± 
as described around (2.13) the explicit matrix form of these terms reads as 



SrpArpB\^ I 
i± 1 2jj |£) ± 



N 

£ 

i=l 



a 



±|2 /rpArpByi 



SjSib ± = diag(|of | 2 ,...,|<7^,0, . ..,0) 



r ±|2 



V -5 AB 



v 2 I N 



(B.1) 



where the Nf x Nf matrix I?j has unit entries on the first diagonal elements and 
is otherwise zero. The expressions on the boundary dD ± of the discs are given up to 
exponentially suppressed terms. The operators on the disc are then given by 

£t£ = -d 2 a + g 2 v 2 l' - (b a d a + c + Ac) , ££t = _^ + g 2 v 2j, _ + c + A ~ c)j (B 2) 

so that they differ only in the last term. The "identity" I 1 = diag(l a , l f if-) has zero 
entries for i, j > N. The respective terms are 



■2i 



Ac 



-Ms 



*t S^l 



(B.3) 



where the off-diagonal terms * ~ T> z T,j vanish exponentially for r — > oo, see (2.12). They 
will drop out in the following analysis. The same applies to c in (B.2), which is diagonal 
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and contains the terms V a A a , (Y< i T A T B Y> i \v ± - ^5 AB ) and (EjEj b ± - v 2 I*j), all of 
which vanish exponentially for r — > oo. Since also the chromo-magnetic field B3 vanishes 
exponentially we note that the two operator products in (B.2) coincide at the boundary 
dD± with only b a vanishing polynomial: 

#D± : CC ] « £ f £ « (-d 2 + 5VI') - Ma , M« ss - — 

r 

where is constant and commutes with <j>o, which can be therefore chosen to be in the 
Cartan subalgebra (x a is the unit- vector) . The term b a therefore vanishes only like ~ - 
and contributes to the surface integral. 



A 1 



e afi x ( 



0,1 



(BA) 



Anomaly. The anomaly contribution in (3.45), with the help of (3.37), can be written as 

^ x) = Tr ^M 2 +% 2M itfC + M* + 4~ Ctf + M* + tf] } ]y=x - (B ' 5) 

As for the overall anomaly, discussed in appendix A, we expand the Green's functions 
in (B.5) for large fj,, see (3.20), keeping the terms that survive the limit /U — > 00 of (B.5). 
The Green's function in question take with (B.2) the form 

(B.6) 



£t£ + M 2 + p? v -d 2 + (g 2 v 2 F + M 2 + fi 2 + p 2 ) - (b a d a + c + Ac) ' 

and analogously for LC) , replacing Ac with Ac. The first term in the expansion is the 
propagator A m (x — y), which is obtained from the D = 2 case of (A.l) by replacing the 
mass (i by the mass-matrix 

m 2 = g 2 v 2 i' + M 2 + fi 2 +p 2 . (B.7) 

Thus the propagator inherits the matrix structure of the diagonal 21 matrix m 2 . However, 
these contributions cancel between the two Green's functions in (B.5). The residual terms 
in the expansion are of the form (A. 5) with Ai — > A m / M and 

0(x +f) = c(x + f L ) + Ac(x + f L ) + Li b a (x + d z « , (B.8) 

and the same for the second Green's function in (B.5) with Ac replaced by Ac. 

On the two-dimensional discs v = and therefore each term in the expansion (A. 5) 
is now regular at y = x, see below (A. 4). Also, the ra'th term in the expansion is now 
proportional to ~ /i _2n+r , r being again the number of derivatives, and thus only the 
contributions (n = l,r = 0,1) and (n = 2,r = 2) survive the limit \i — > 00 of (B.5). 
However, the contributions with the maximum number r = n of derivatives are the same 
for both Green's functions in (B.5) and thus cancel each other. 



21 Note that the entries of M, see (3.32), live in the Cartan subalgebra so that the fundamental and 
adjoint entry of M 2 are diagonal, see appendix C. 
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Therefore one is left with the single contribution (n = l,r = 0) to (B.5): 



a(x) = Tr { 



M 2 



M 2 + M 2 



2 Mil' 2 A m/M (-z) [Ac - Ac] (x + *) A m/M (z) } . . . 



Tr 



7r I .M 2 + M 2 5 2 v 2 /' + .M 2 + n 2 + p 2 



(B.9) 



In the second step we evaluated the difference [Ac — Ac] at x and used that at the discs 
D± the field B commutes with Ai, and thus with the propagators A m / M , see below (2.13). 
In addition, we note that A^ ft (0) = 4^^, see (A. 3) with v = 0. 

The actual local anomaly contribution (3.45) is the integration over p of a(x). However, 
it commutes with the limit [i — > 00 which puts the second factor in the trace to one. This 
leaves for the p-integration, 



M 



M 



dp 

2tt {M + M 2 +p 2 ) ~ 2VM 2 + M 2 



(B.10) 



where we note that the integral is well defined. The possible zero eigen- values of M 2 drop 
out because of the factor Ai in the numerator. Therefore even for M — > 0, a limit that is 
eventually of interest, the integrand has no poles. 

We thus obtain for the disc anomaly contribution to the index function (3.27), (3.45), 



A 



disc 



t7 M 



M 



VM 2 + M 2 



(B.ll) 



which is in fact M -dependent. 



Surface term. Contrary to the anomaly, the surface term is evaluated at the boundary 
dD± of the discs. Using (3.37) the current in (3.45) can be written as 



q a (x) = ±Tr {TT a [n 



1 



-un 



y=x , 



(B.12) 



where the operators involved are defined in (3.36) and the last term indicates the sub- 
traction of the same expression with M — > fi. In the following we will first look at the 
expressions for M. Up to exponentially suppressed terms LC) and C) C are identical at 
the boundary of the discs (B.4) and therefore the inverse operator in (B.12) considerably 
simplifies close to the boundary, 



-un 



i 4 ® 



-d 2 + m 2 M 



b a d a ' 



(B.13) 



where we have again introduced a mass matrix m 2 ^, which is the same as in (B.7) but \i 
is replaced with M, as indicated by the index. At the boundary b a vanishes like ~ r , 
see (B.4) and it is sufficient to know the Green's function (B.13) up to orders 0(r~ 2 ) to 
perform the integral over dD±. 
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The expansion of (B.13) is of the form (3.20) with A M replaced by A mju and the 
insertion operators are O(z) = b a (z) d z <x. Changing the integration variables as Zi — > x + z% 
the insertion becomes b a (x) d z a + 0(r~ 2 ) and the Green's function (B.13) is approximated 
by 

A mM (x -y) + b a (x) d xa A 2 (x -y) + 0(r~ 2 ) , (B.14) 



b a d a 



where, as in the case of the disc anomaly, we have used that the gauge fields and thus b a 
(B.3) commute with the mass-matrix m 2 M and the propagator A TOm , see below (B.9). 

We note that the product YY a in (B.12) is block off-diagonal, see (3.38), (3.42), whereas 
is block diagonal. Consequently, only the off-diagonal terms of % contribute in the 



-HH 



trace (B.12), which are independent of M. We thus have 

q a (x) = \ Tr { rr Q ^ 1 4 ® ( [A mM (x - y) - A m ^(x - y)] 



+ b a (x) d xa 



A 2 



( x -y)-Al^x-y)])}\ y=x + 0(r- 2 ) 



(B.15) 



The expressions in the two square brackets are regulated and well defined at y = x. We 
also note that T ~ jvP+Af 5 ' w ^^ cn ^ s an important convergence factor for the p-integration 
(3.45) and as in the case for the anomaly projects out possible poles for M = 0, which 
appear here also in the propagator because I' in the mass- matrices of the form (B.7) does 
not have full rank, see below (B.2). We can therefore safely set y = x and the resulting 
current is then, 



q a {x) 



Tr { TT a K(x) U ® log 
1 



m 



M 



tnf 



} 



f — Tr{rr Q Tpl 4 ®b (x) log 

1D7T l 



m 



mi 



}+0(r~ 2 ) 



(B.16) 



For the first term we used that the derivative part of % (3.38) vanishes by symmetric 
integration. Similarly, single derivatives vanish for the second term and <9 a <9^A 2 eg = 



2 



5 n nd 2 A 2 



Seemingly the limit \i — > oo of (B.16) does not exist, but this is the last step to carry 
out in (3.45) and we note that after the p-integration the limit does exist. However, before 
doing so we use that I4 <g> b a = {Tp, K} and that the matrix log [m^/m 2 ] commutes with 
T a such that with the relations (3.40) one can show that the two terms in (B.16) in fact 
cancel. Thus q a (x) = and the disc surface term does not contribute, as stated in the 
main text. 



C. Cartan-Weyl basis for U(N) 

We give some basic definitions for the Lie-algebra of U(N) and thereby introduce the 

conventions that are used in the main text. Throughout we use the notation e^i jy for 

the canonical basis of the vector space M. N , i.e. with unit entry at position m and otherwise 
zero (in components e l m = 5^) 
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In addition to the rank = N — 1 mutually commuting generators of su(iV) one has 
the unit matrix (2.4) in the Cartan subalgebra f) C u(N). The absence of the traceless- 
condition allows the choice of a rather convenient basis for f). Together with the ladder 
operators these form a basis of u(N): 

(C.l) 

The N Cartan generators Hi have the non-vanishing entry at the i'th position on the diag- 
onal, whereas the N(N — 1) ladder operators E amn have the only non-vanishing entry at 
(row, column) = (m, n) with m ^ n. The ladder operators are labeled by the ^-component 
roots a mn = (a^n)) which are defined by the second of the following commutation rela- 
tions: 22 



H 



1 

72 



[Hi,Hj] = 0, [Hi,E a ] = a l E a , [E a , E„ a ] = a 1 H t . (C.2) 
The N(N — 1) roots are vectors G M. N and are of the form 

The set of roots S := {a mn } can be divided into positive and negative roots, E = S + U S~ 
by the convention that a root is positive if the first non-vanishing entry is positive and 
otherwise negative. The positive roots are therefore given by S + = {a m<n }. The negative 
roots are just the negative thereof. Also the positive roots are not linearly independent. A 
set of linear independent roots that generates all roots by linear combinations with integer 
coefficients of a common sign is called simple roots: 

- n— 1 

(3 a = —={e a - e a+1 ) , a = l,...,N-l S ± B a mn = ± ^ Pa , (C4) 

* a=m 

thus positive/negative roots have positive/negative integer coefficients. 

From the definition of the ladder operators and the roots follows that E& = E^ a . By 
defining the hermitian generators 

Tq + ' = -j=(E a + E-a) , Ta ' = —jJyEa — E^ a ) , (C5) 

a hermitian basis of generators with the normalization as given in (2.4) is then given by 
{Hi,T& ,Ta ' | a € These relations are completely standard and are the same 

as for 5u(N) except for the fact that the Cartan subalgebra has now one more element. 
Basically the only choice of conventions that matters is the normalization (C.3). 



22 Here and in the following we often omit the labels m, n for the different roots. A root symbol a, (3 etc. 
stands then for a particular root, i.e. particular values for m, n. 
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It will be convenient to have an explicit matrix realization of the Cartan generators in 
the adjoint representation. In the hermitian basis introduced here it is given by 



trad 







iA k 
-iA k 



(C.6) 



where A k is the diagonal matrix of the fc'th component of all positive roots, i.e. A k 
diag(a^ 2 i • • • > a %-\ at)- This implies the normalization 



Tr ad {HiHj} = 2 J2 a * « j = Yj ai aj = N ^ ■ 



(C.7) 



Note, in the non- hermitian basis (C.l) the adjoint matrix of Hi is diagonal. 
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